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Complex spatiotemporal structures develop during the process 
of aging glasses after cooling and of rejuvenating glasses upon 
heating. The key to understanding these structures is the inter- 
play between the activated reconfiguration events which generate 
mobility and the transport of mobility. These effects are both ac- 
counted for by combining the random first order transition theory 
of activated events with mode coupling theory in an inhomoge- 
neous setting. The predicted modifications by mobility transport 
of the time course of the aging regime are modest. In contrast, 
the rejuvenation process is strongly affected through the propa- 
gation of fronts of enhanced mobility originating from the initial 
reconfiguration events. The structures in a rejuvenating glass re- 
semble flames. An analysis along the lines of combustion theory 
provides an estimate of the front propagation speed. Heteroge- 
neous rejuvenation naturally should occur for glasses with free 
surfaces. The analogy with combustion also provides a new way 
of looking at the uptake of diluents by glasses described by case 
II and super case II diffusion. 

Introduction 

It is an excellent, but still imperfect, approximation to think of a 
structural glass as being a frozen snapshot of the liquid state. The 
macroscopic properties of a glass depend on its preparation history in 
contrast to those of an equilibrated liquid. Likewise, the structure of a 
glass is history dependent and continues to change even as it is being 
studied. The picture of the dynamics of an equilibrium supercooled 
liquid that emerges from the random first order transition (RFOT) 
theory of the glass transitionfl] can be described, with some poetic 
license, as that of a fluctuating mosaic of structures that locally are 
chosen from a diverse set of minimum energy patterns which recon- 
figure through activated events. In the nonequilibrium aging glass, 
this mosaic does not just fluctuate but continues to evolve, eventually 
to the dynamic mosaic of an equilibrated liquid sample at a lower 
temperature[2]. An aging glass can also be heated up, "rejuvenating" 
the sample, again eventually reaching a new dynamic equilibrium 1 . In 
this paper, the spatiotemporal structure of the glassy mosaic in the in- 
termediate states of development of aging or rejuvenating is explored 
within the framework of the RFOT theory. 

A spatially ultra-local description of an aging glass has already 
been provided by Lubchenko and Wolynes[2]. As in equilibrium, 
activated transitions of small regions in the nonequilibrium glass are 
driven by the extensive configurational entropy of local structural pat- 
terns. Such transitions (sometimes called "entropic droplets") are 
retarded by the stability of the local pattern, which varies spatially, 
giving a distribution of relaxation times. The range of existing lo- 
cal stabilization energies in a small region may be described by a 
Lagrange multiplier, the Active temperature which can be thought of 
as fluctuating spatially [3]. The least stable local regions typically re- 
configure first. Upon cooling, the mosaic, therefore, begins to contain 
patches of greater stability than it had before. In this way, the aging 
glass not only on the average becomes more stable, but also, for a time, 
is even more inhomogeneous in its energy distribution than an equi- 
librated sample since it was the fast patches that became slow first. 
The average behavior leads to non-linear relaxation much like that 
used in the Nayaranaswamy-Moynihan-Tool phenomenology which 
is based on a uniform Active temperature[4, 5, 6]. The results of the 
microscopic RFOT theory for the nonlinearity parameter in the NMT 



approach are in good agreement with experiment. In addition, how- 
ever the RFOT theory also predicts a qualitative difference from NMT 
phenomenology. RFOT theory suggests the patchier mosaic found at 
intermediate times should lead to additional "ultra-slow" relaxations. 
These anomalous relaxation processes have been observed in some 
experiments [7]. 

The local picture of aging does not completely account for the dy- 
namical coupling between re-configuring regions: the reconfiguration 
of a given region changes the pinning forces acting on its neighbor- 
ing regions, allowing their reconfiguration rates, in turn, to change. 
This effect, may be termed facilitation, as in the popular kinetically 
constrained models of glasses[8, 9], where it is typically the only ef- 
fect being modeled. Facilitation was noted as an aspect of aging by 
Lubchenko and Wolynes but was not completely analyzed by them. 
Bhattacharya et al. (BBW) have shown that the facilitation of acti- 
vated processes in equilibrium liquids can be captured mathematically 
by combining mode coupling theory (MCT) with the entropic droplets 
that describe activated transitions inRFOT theory[10, 11]. Weexploit 
this insight to describe the nonequilibrium spatiotemporal aspects of 
glasses by using the combination of MCT with the RFOT mosaic to 
motivate a continuum theory for the spatiotemporal coupling between 
activated events in a glass or supercooled liquid. An equation for what 
may be called a mobility field[ 1 2 ] , follows from inhomogeneous mode 
coupling theory within this BBW-based framework. Unlike existing 
treatments of inhomogeneous mode coupling theory, the equation has 
a spatially varying source term. This source arises from the activated 
events, whose dynamics within RFOT depend on the local energy 
or Active temperature. The local energy, in turn relaxes to local- 
equilibrium at a rate that depends on the local mobility at the same 
location. The resulting coupled equations for a mobility held and for 
a local energy or Active temperature resemble the nonlinear diffusion 
equations encountered in the theory of combustion[13, 14]. In com- 
bustion, chemical kinetics depends strongly on temperature which 
is increased through heating from the reaction events themselves, but 
that also is transported by conduction. For the aging glass, the coupled 
nonlinearities will be shown to have a quantitative but not qualitative 
effect on the spatiotemporal evolution of the mosaic. On the other 
hand, the situation is quite different for the rejuvenating glass. The 
rejuvenating glass, like a highly combustible mixture is unstable and 
once reconfiguration events are nucleated they can propagate through 
a sequence of additional reconAguration events, as do exothermic re- 
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1 The term "rejuvenation" is used here as it is in the spin glass community. The polymer com- 
munity uses the term to describe the stress induced reversal of aging. This latter process is 
distinct from but may be related to the present heating protocol. 
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action events in a flame. These flame structures in the glass speed 
up the rejuvenation process enormously. While specific quantitative 
evidence for such structures has not yet come to light in a homoge- 
neously rejuvenating glass, there are hints from imaging methods like 
those pioneered by Israeloff[15]. Ediger has observed front propaga- 
tion that initiates at the surface of vapor deposited ultrastable glass 
upon heating it[16, 17]. The speed of the observed heterogeneous 
rejuvenation front seems to be consistent with the predictions of the 
present theory. 

The plan of the paper is as follows. A route to the equations for 
the mobility field driven by activated events is described. An esti- 
mate for the effect of facilitation in the ordinary aging situation is ob- 
tained. Front propagation in the rejuvenating glass is then treated us- 
ing some of the simpler approximations from combustion theory [14]. 
This simple treatment allows an estimate of the rejuvenation time for 
a homogeneous sample without surfaces. The origin of heteroge- 
neous rejuvenation within RFOT is discussed. Finally the coupling 
of these spatiotemporal processes to plasticizer diffusion suggests an 
explanation within the RFOT theory framework for the mysterious ac- 
celerating penetration of diluents into polymers known as Supercase 
II diffusion! 18]. 



Mobility Transport and Mobility Generation in Glasses 

While mode coupling theories are often formulated in a momentum 
space representation appropriate for a uniform system, the physics 
underlying the equations is as valid for inhomogeneous systems as it 
is for homogeneous ones. Long ago the nonlinear coupling of hydro- 
dynamic modes in a fluctuating simple fluid with surfaces was shown 
to modify the hydrodynamic boundary conditions [19]. Coupling of 
structural density fluctuation modes in an inhomogeneous situation 
also has been shown to describe a growing length scale as the transi- 
tion to nonergodic behavior is approached[20, 21]. In structural glass 
dynamics, mode coupling theory relates the memory kernel for lo- 
cal structural relaxation to the behavior of density-density correlation 
functions. These correlation functions in turn depend on the same 
memory kernel but, in an inhomogeneous system at locations dis- 
placed in space and time. The memory kernel is frequency dependent, 
for time translation invariant systems. This frequency dependence is 
interesting and complicated, leading to behavior resembling (3 relax- 
ation. In the present paper we will, however, simplify the analysis 
by treating the memory kernel for the explicitly time dependent ag- 
ing or rejuvenating systems, as a time and space varying, but locally 
frequency independent, rate, /x(r, t). This rate can be called a mobil- 
ity field. Biroli et al., generalized ideal mode coupling theory to the 
inhomogeneous situation by expanding the self-consistent equations 
in the degree of inhomogeneity and in gradients to obtain a contin- 
uum description[21]. Carrying out this procedure for the memory 
kernel of mode coupling equation containing an activated event term, 
as proposed by BBW, yields an equation of the form 

Here p is the uniform solution of the homogeneous MCT equations 
with activated transitions. In the strongly glassy regime MCT effects 
are modest. Therefore p will depend on the local temperature and 
Active temperature, in a way very close to the existing Lubchenko- 
Wolynes (LW) theory for uniform systems that accounts for activated 
transitions alone. 

By taking the static limit we can see that the coefficient of the 
linear Laplacian term contains a length £ which is the correlation 
length of the 4-time correlation function consistent with the Biroli et 
al. analysis. The gradient squared term arises because of the nonlin- 
ear relation of MCT closure between memory kernel and correlation 
functions. The coefficients £ 2 and c, depend on the details of the mi- 



croscopic mode coupling closures employed. For simplicity we will 
choose c so that the locally linearized equation can be written as a 
mobility flow equation with a source 

^ = - v 'Jm-m(m-m) [2] 
ifi = - mC 2 Vm 

In the nonequilibrium glass we must also take the local stabilization 
as a variable-equivalently we may say there is a spatially varying 
Active temperature that satisfies an ultra-local relaxation law without 
conduction when the finite spatial structure of the activated events is 
neglected: 

= - T) [3] 

The actual temperature equilibrates by vibrational thermal conduc- 
tion rapidly and will be treated as uniform. In the absence of mobility 
transport through a flux term, these equations are equivalent to the 
Lubchenko-Wolynes version of the NMT formalism. 

To describe the fluctuating Active temperatures and energies 
known to occur in the equilibrium liquid these equations for the mean 
behavior must be supplemented by random force terms that are lo- 
cal. Owing to these random forces the mosaic will fluctuate before a 
quench or heating process is initiated, then evolve from these initial 
fluctuations and then proceed to fluctuate again. The intensity of the 
fluctuations in the mobility transport equation reflects the shot noise 
of creating entropic droplets while the fluctuating forces in the Active 
temperature equations ensure there is a trend to proper Auctuating 
thermal equilibrium. For the situations considered here knowledge of 
the intensities, but not of the precise statistics, of these Auctuations is 
sufAcient. 

While p in pure mode coupling theory undergoes a divergence at 
the nonergodicity transition, due to the activated events that cut off 
the transition in RFOT theory, p instead becomes a strong, super Ar- 
rhenius function of the Active temperature. The strong dependence of 
p on Active temperature is the key to the analogy to combustion for 
rejuvenating glasses where reaction rates also depend exponentially 
on the local temperature. The Lubchenko-Wolynes RFOT analysis of 
activated events in a nonequilibrium system suggests that the tradi- 
tional NMT formulation of p using locally two distinct Arrhenius laws 




Fig. 1. The mosaic pattern of mobilities soon after a quench from the tempera- 
ture To to a lower temperature T, as the sample begins to age. A few regions 
of size £o have equilibrated to low fictive temperature structures with 7> to near 
to T, while others maintain their higher initial fictive temperature. The initially 
equilibrated regions now have a low mobility /ji ow , contributing to the "ultraslow" 
relaxations. Mobility transport to these regions also slows approach to equilibrium 
of the neighboring cells but the effect is modest. 
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for Active and ambient temperature[2] should be reasonably adequate 
if the quenches are not too extreme: 

fi(T F , T) = exp | ~ K T J j [4] 

When T — Tf this reduces to the Adam Gibbs form, according to 
RFOT theory. The more complete LW RFOT expression containing 
the configurational entropy explicitly could equally well be used. In 
RFOT owing to the configurational entropy dependence the nonlin- 
earity parameter x is predicted to be a function of the configurational 
heat capacity AC p. This result of RFOT theory agrees well with 
experiment[2]. The BBW-MCT for a uniform system should give 
only mild modifications to the local LW-RFOT relation due to the av- 
erage effects of facilitation while the spatiotemporal structures con- 
sidered here from the inhomogeneous version of the theory can have 
qualitative effects on the global process of equilibration. 



Aging Glasses 

How is the spatiotemporal pattern in the aging glass after rapid 
quenching modified by the mobility transport that was not studied ex- 
plicitly by LW? Consider a rapid quench to an ambient temperature 
T from an equilibrated sample at a temperature To which is then at 
the average initial Active temperature TJ? = To but that also has local 
Tf fluctuations given by the configurational heat capacity. The ini- 
tial dynamics does not involve mobility transport but only activated 
transitions as already described by LW. This initial step generates 
droplets of radius £o> that are nearly equilibrated. The radius £o is 
predicted by LW to typically be about 5 particle spacings, rn, near 
the laboratory glass transition, T g . These initial transitions occur at 
random, starting in those regions where Tf had fluctuated to a higher 
than average value. After a reconfiguration event, each of these re- 
gions has come to equilibrium at a Active temperature near to the 
ambient temperature T. As shown in figure the mosaic, after some 
initial transitions has occurred, is more inhomogeneous than a sample 
characterized by a Active temperature with only Gaussian fluctuations 
would be. The number of the reconflgured regions per unit volume 
during this initiation phase is thus (ro/£o) 3 (l — <l>(t, Tf, T)) where 
4> is the relaxation function of the glass. Neglecting the possible role 
of an early f3 process, the relaxation function for either a supercooled 
liquid or glass within RFOT theory is tolerably well approximated as 

a stretched exponential e - ^ '' where fio = fi{T F = To, T) 

Each of the regions that have reconfigured now becomes ultra- 
slow, since each one typically has a greater stability than it had be- 
fore. This is the origin of the anomalous relaxation component found 
by Miller and MacPhail[7]. Mobility transport now will affect the 
surroundings of each of the initially reconfigured regions. These sur- 
rounding regions relatively speaking, reconfigure more rapidly than 
the newly equilibrated regions, being at the higher initial Active tem- 
perature. The initial reconfigured regions can thus be considered a 
static influence on the relaxation of their surroundings. Nevertheless 
if T is considerably lower than To, then the mobility of the remain- 
ing regions to be equilibrated is not just £l(Tf,T) but is additionally 
slowed by mobility transport to the ultraslow inclusions in the mosaic 
which have mobility jUi ow . As shown in figure each of the low Active 
temperature regions affects a region around it, lowering the mobility 
following the quasi-static law [22] 

-C 2 VV+(M-M) = [5] 

where p, is very small within a previously reconflgured region but 
large elsewhere, fi is decreased significantly only within a length £ of 
a drop. It follows that even in the perturbed surroundings of an ultra- 
slow cell the rate typically is changed through mobility transport by a 
modest factor of roughly two until the reconfigured regions are very 




Fig. 2. The mosaic pattern of mobilities for a sample originally equilibrated at 
the temperature To which has just been raised to the temperature T, beginning 
its rejuvenation. The initially nucleated regions of high mobility catalyze the rear- 
rangement of their low mobility neighbors leading to a radially propagating front 
moving at a velocity v. 



close to merging. The regions of influence of the initially reconfigured 
mosaic cells will overlap when 

(2£) 3 C - 3 (W) ^ 1^0^7/8 [6] 
The stretched exponential relaxation implies 

^overlap ~ 77^0 t^l 
O 

Rejuvenating Glasses 

The initial stage of rejuvenation of a glass proceeds upon heating 
very much as does aging after a quench but with two important dif- 
ferences. Since the ambient temperature is increased, rather than de- 
creased, the first reconfiguration events occur faster than they would 
in the original equilibrated sample (not slower as in aging). These 
events still do not happen quite as fast as they would occur in a sam- 
ple already equilibrated to the higher ambient temperature since the 
sample is initially more stable than an equilibrated high temperature 
sample would be. Since the reconfigured regions are now at high Ac- 
tive temperature internally their motions become faster. Indeed their 
mobilities are exponentially larger (not smaller!) than the mobility 
of their surroundings as shown in Agure . Mobility transport now 
allows the immediately surrounding material around an initially re- 
configured region to begin to change more rapidly than it would on 
its own. The mobility increase in neighboring regions is autocatalytic 
and therefore a front of higher Active temperature should emanate ra- 
dially from each initially rejuvenated center. These growing zones of 
infiuence will quickly overlap and the glass, as a whole will rapidly be 
equilibrated. This situation is shown in Agure . Of course, this simple 
scenario is most appropriate when the heating jump is large, which 
is the case we will explicitly analyze. The time when overlap of the 
rejuvenating regions is reached will be tn, termed the rejuvenation 
time. 

Let us assume the front can propagate stably. The front moving 
at a velocity v sweeps out at a volume (wt) 3 in a time t. Overlap will 
be reached when 

1~ (vtRfZo 3 (l-0(t R ,T F \T)) [8] 
Again approximating by a stretched exponential we obtain 
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Fig. 3. A sketch of the mobility front in a comoving frame. The mobility rapidly rises from /ii ow , before the region has equilibrated, to /ihigh ~ ft(T) after a region 
has rejuvenated. The inverse fictive temperature Tp 1 falls while /i rises. This figure parallels the structure of a flame front, where, however, the temperature in the 
flame is the analog of mobility here and fuel concentration is the analog of Tp 1 . Mobility generation occurs where there is both a significant disequilibrium in mobility 
while the mobility is still rather high. 
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[9] 



This implies tR is essentially a weighted geometric mean of the time 
it would take the front to cross a nucleated drop and the original re- 
laxation time for a system starting at the initial fictive temperature Tf 
but ambient temperature T: 



3/(3+/3) 



{r(Ti?,T)} 



0/(3+13) 



[10] 



The analogy with combustion allows an estimate of the front velocity. 



Rejuvenation Front Propagation and the Combustion Anal- 
ogy 

Consider the coupled mobility/fictive temperature transport equa- 
tions for a planar front in which /i and Tf vary only along a single 
dimension, x. This situation is shown in figure 3. As in the theory 
of a flame front[14], the coupled equations can be written in a frame 
comoving at velocity v in terms of one variable z 



V ~dz ~ ^ 



■fi) = 



[11] 



T) — 



where jl — JX{Tf,T). In the original frame of reference both Tf and 
/i are functions of z = x + vt. /i and Tf are graphed schematically 
in figure 3. We consider the case where the mobility at the rejuvena- 
tion temperature is much larger than the original fictive temperature. 
According to these equations mobility is generated to a significant de- 
gree only in a thin generation zone localized at the front near z = 0. 
Outside this zone either the fictive temperature has already had time to 
equilibrate so /i = ft = /ihigh or else fj, itself is small because the sam- 
ple is too slow to have equilibrated at all yet so the mobility generation 
rate is low. The latter condition (no equilibration in the unperturbed 
region) is familiar in combustion theory where, in simpler treatments, 



it is usually assumed that the unreacted mixture is completely stable. 
This assertion is not strictly true since the mixture would eventually 
react but on an exponentially long time scale in the absence of auto- 
catalytic heating by a pre-existing flame. In the current context this 
infamous "cold boundary difficulty" is not a problem under the con- 
dition that rejuvenation jump is sufficiently large. Flame propagation 
and front propagation for a random first order transition correspond to 
a problem of "intermediate asymptotics" in the language of Barenblatt 
andZel'dovich[23], 

Since z does not explicitly enter into equation 11, it is conve- 
nient to solve the equations by introducing explicitly the mobility flux 
j — —^fld^/dz and to re-write equation 1 1 as an equation for j as 
a parametric function of /i itself, leaving the z dependence implicit. 
This gives 



3^+ v 3 - Ma» - = 



[12] 



Outside the mobility generation zone the first two terms involving 
mobility transport balance giving 



vfi + constant 

V{H - //low) 



[13] 



As approaches /ihigh, it becomes possible for the structural rear- 
rangements to take place rapidly enough to lead to mobility gener- 
ation. In the mobility generation zone, a different approximation to 
equation 12 is valid in which the first and third terms of equation 12 
are balanced rather than the first two. The balance in the generation 
zone gives an alternate formula for the mobility flux in the comoving 
frame: 



Idf _ 2 2 
Integrating this equation at the high mobility side yields 

11/2 



/-Mhigh 

2£ 2 / 



= j (Mhigh) 



[14] 



[15] 



Matching this with the crossover to the mobility transport region gives 
an expression for the front velocity 
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/^high A*low 



< J Clow 



1/2 



[16] 

To evaluate this expression exactly the Active temperature profile 
needs to be known i.e. how quickly 7> varies in the reaction zone. 
The Tf profile reflects how far the system is dragged out of local equi- 
librium through mobility transport versus local relaxation to the new 
local thermal equilibrium. This requires the simultaneous solution of 
the Tf and fi equations and must be done numerically. A fairly accu- 
rate estimate (an over-estimate in all likelihood) for the propagation 
speed can be found, however. 

This estimate for the velocity follows from the fact that the dis- 
equilibrium of mobility fj, — p,(T F ) must be smaller than, but can 
be of the order of, the overall change in mobility across the front 
Mhigh — Miow Introducing this overestimate of the disequilibrium in 
the mobility generation zone into equation 16 gives then 



(/Xhigh — Mlow) 1/2 



/-Mhigh 

2f / dw 2 {T F ) 

J Maw 



1/2 



[17] 



In this expression the behavior near the high mobility end would dom- 
inate. At the high mobility end p,(T F ) ~ fi. This relation can also not 
be too far wrong at low fi where again fi must settle down to /n(Ti ow ). 
Using the approximation here that fi — jl uniformly in the integral in 
equation 17 yields an estimate valid for /Lthigh S> /J,\ow- 

v = x/^/ihigh [18] 

When the overall mobility change is large, the width of the mobility 
generating zone v/fi is thus predicted to be of the order £, not too 
different from the size of an entropic droplet. This scale itself is only 
modestly larger than the fundamental molecular size. This suggests, 
as in the theory of shock waves where the shock width is of the order 
of a mean free path, a more complete microscopic treatment, going 
beyond the continuum treatment may be necessary for quantitative ac- 
curacy. The microscopic inhomogeneous BBW-MCT can be solved 
but it appears to be a daunting numerical task. Indeed, making the 
approximation of a molecularly sharp interface would not be entirely 
out of the question for analysis. Such a sharp interface would give a 
rate of roughly the same magnitude. 

We have so far neglected fluctuations entirely in treating front 
propagation. Such fluctuations may give rise to instabilities and more 
complex front structures. Nevertheless it is the faster reconfiguring 
regions that should matter the most. To a first approximation fluc- 
tuations should lead to a velocity which averages the mobility of the 
rejuvenated sample rather than averaging the reconfiguration time. 

Combining the estimate from equation 18 with the earlier equa- 
tion for the rejuvenation time one obtains 



t R 



14) 



3/(3+0) 



[f(T|'\T)] 3/(3+/3) [r(T F \T)] 



0/(3+0) 



[19] 

Here f is the harmonic mean relaxation time at the final Active tem- 
perature Tp n ~ T. Since £ ~ £o the prefactor is of order unity. 
We thus see the rejuvenation tu is a weighted geometric mean of the 
initial and final relaxation times. We have assumed in the analysis of 
the front that the mobile phase has equilibrated after the passage is 
complete. It is therefore useful to notice that this relation confirms 
that the assumption of stationarity is indeed justified since the initial 
relaxation time is indeed much longer than the final one. 

More weakly rejuvenated samples (i.e., T greater but close to T'p) 
cannot be analyzed with the simple constant velocity formula for front 
propagation. On the other hand a linearized analysis of the coupled 
mobility and Active temperature equations should be adequate to treat 
such cases. 



Heterogeneous Rejuvenation 

While aging of a macroscopic sample should be dominated by pro- 
cesses in the bulk, the rejuvenation of a glass should be faster at the 
surface than it is in the bulk because the activated processes that orig- 
inate mobility generation are faster at the surface. Stevenson and 
Wolynes have shown that RFOT theory implies a reduction of the 
reconAgurational activation free energy by a factor of two at a free 
surface in an equilibrated liquid and the same factor should apply to 
glasses in the aging regime[22]. Mobility generation will hence pro- 
ceed much more rapidly at such a free surface and a rejuvenating front 
will preferentially start at the surface and propagate into the bulk. Such 
heterogeneous rejuvenation was observed by Ediger's group starting 
with an ultrastable glass[17]. They describe the observed speed of 
front propagation as having a temperature dependence paralleling the 
rate of molecular self diffusion. The present arguments are consistent 
with their observations since the self-diffusion constant in a heteroge- 
neous system averages reconAguration rates not times, just as the front 
velocity here is related to average rates of structural rearrangement. 

The stability of front propagation, a central focus of so much com- 
bustion science, may not be too critical for homogeneous rejuvenation 
since little time elapses before the fronts merge. The stability analysis 
is expected to be much more relevant to the heterogeneous rejuvena- 
tion starting from the surface. Since the effective temperature is not 
directly transported in the present model, the situation might seem 
to resemble combustion of condensed media more than the combus- 
tion of gasses. On this basis, front propagation would be expected to 
be stable but clearly a detailed mathematical stability analysis would 
be required to establish this, especially when fluctuations in Active 
temperature are explicitly included. 

Front propagation from surface to bulk is a well-studied feature 
of diluent penetration into glassy polymers, where it is referred to as 
case II diffusion [18] in which uptake is linear in time, rather than 
the square root expected for simple diffusion. A large number of 
theories of case II diffusion exist that have many similarities to the 
present study in purely mathematical terms[24, 25]. There is a major 
difference, however, between those treatments and the present one. 
Even without a diluent, that can diffuse, the aging and rejuvenation 
phenomena studied here by themselves lead to a diffusion equation 
for the mobility. If at the same time a diluent diffuses into the glass, a 
double diffusion problem appears. If the diluent molecule is smaller 
than the molecules in the bulk of the glass, it can diffuse quite rapidly, 
thus there will be considerable decoupling of the diluent diffusion 
time scale from the glassy a relaxation time which determines the 
diffusion of mobility. In combustion the ratio of a reactant's diffu- 
sion coefAcient to the kinematic thermal conductivity is known as the 
Lewis number[14, 13, 26]. The corresponding ratio of diffusion coef- 
Acient to mobility diffusion coefAcient L M will be important for type 
II diluent penetration. 

A large Lewis number in combustion leads to unstable, acceler- 
ating flame front propagation often leading to eventual thermal deto- 
nation. In analogy, if the diluent diffusion is strongly decoupled from 
glassy relaxation i.e. if L M is large, the penetration front will acceler- 
ate, leading to more rapid than linear uptake of the diluent. Although 
a more complete mathematical analysis than is done here is needed, 
the analogy with combustion evident from the present framework thus 
provides a natural explanation of accelerating uptake of penetrating 
diluents. Such accelerating uptake has been observed and is called 
supercase II diffusion. 



Conclusion 

The combination of the generation of mobility by activated reconAgu- 
ration events and mobility transport leads to complex spatiotemporal 
structures in glasses. Aging and rejuvenating are seen to be, in no way 
reciprocal processes mechanistically[27], although they are described 
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by a common set of equations. More complex thermal histories can 
give rise to more complex spatiotemporal patterns than those described 
here for simple heating and cooling experiments. It is hard to believe 
these patterns are not an essential aspect of many heat treatment pro- 
tocols used in technology, but clearly, as in the study of combustion, 
the range of phenomena can be quite rich and difficult to analyze 
mathematically, even when the best engineering practice is clear em- 
pirically. Only the simplest cases have been discussed in this paper. 
In any event, the present analysis encourages us to contemplate the 
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